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Abstract. We present analyses of the two-point corre- 
lation properties of the ESO Slice Project (ESP) galaxy 
redshift survey, both in redshift and real space. From the 
redshift-space correlation function £(s) we are able to 
trace positive clustering out to separations as large as 
50 h _1 Mpc, after which £(s) smoothly breaks down, cross- 
ing the zero value between 60 and 80 h _1 Mpc. This is best 
seen from the whole magnitude-limited redshift catalogue, 
using the J3 minimum-variance weighting estimator. £(s) 
is reasonably well described by a shallow power law with 
7 ~ 1.5 between 3 and 50 h" 1 Mpc, while on smaller scales 
(0.2 - 2h~ 1 Mpc) it has a shallower slope (7 ~ 1). This 
flattening is shown to be mostly due to the redshift-space 
damping produced by virialized structures, and is less evi- 
dent when volume-limited samples of the survey are anal- 
ysed. 

We examine the full effect of redshift-space distortions 
by computing the two-dimensional correlation function 
£(r p ,Tr) , from which we project out the real-space £(r) 
below 10 h -1 Mpc. This function is well described by a 
power-law model (r/r )~ 7 , with r Q — 4.15t ; 21 /i 1 Mpc 



and 7 = 1.67tg;{g. 

Comparison to other redshift surveys shows a consis- 
tent picture in which galaxy clustering remains positive 
out to separations of 50h _1 Mpc or larger, in substantial 
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agreement with the results obtained from angular sur- 
veys like the APM and EDSGC. Also the shape of the 
two-point correlation function is remarkably unanimous 
among these data sets, in all cases requiring more power 
above 5h _1 Mpc (a 'shoulder'), than a simple extrapola- 
tion of the canonical £(r) =(r/5) -1,8 . 

The analysis of £(s) for volume-limited subsamples 
with different luminosity shows evidence of luminosity seg- 
regation only for the most luminous sample with Mt,j < 
—20.5. For these galaxies, the amplitude of clustering is 
on all scales > 4h _1 Mpc about a factor of 2 above that 
of all other subsamples containing less luminous galax- 
ies. When redshift-space distortions are removed through 
projection of £(r p , tt) , however, a weak dependence on lu- 
minosity is seen at small separations also at fainter magni- 
tudes, resulting in a growth of r Q from 3.45toJo h _1 Mpc 
to 5.15lg' 3 4 h _1 Mpc, when the limiting absolute magni- 
tude of the sample changes from M = —18.5 to M = —20. 
This effect is masked in redshift space, as the mean pair- 
wise velocity dispersion experiences a parallel increase, ba- 
sically erasing the effect of the clustering growth on £(s) 



keywords - Cosmology - Large-Scale Structure 



2 



L. Guzzo et al.: Clustering in the ESP Galaxy Redshift Survey 



1. Introduction 

The spatial two-point correlation function, £(r) , is prob- 
ably the most classical statistic used in cosmology for 
clustering analyses. Since its early applications (e.g. Pee- 
bles 1980), obtaining reliable estimates of on large 
(> 5 — 10h -1 Mpc) scales became one of the main sta- 
tistical motivations for enlarging the available 3D samples 
through new, wider and deeper redshift surveys. The prin- 
cipal reason for this is that on scales sufficiently large for 
the fluctuations to be still in the linear clustering regime, 
the shape of £ (r) is expected to be preserved during grav- 
itational clustering growth. Comparison of observations 
with models is therefore simpler, because in this case the 
theoretical description of clustering does not require the 
full non-linear gravitational modelling, which is on the 
contrary necessary at small separations. If we are able to 
measure accurately £(r) [or its Fourier dual, the power 
spectrum P(k)], on large enough scales, we shall have a 
measure of the distribution of the initial fluctuation ampli- 
tudes. In addition, if the initial density held was described 
by a Gaussian statistics, this will be all the statistical in- 
formation we need to completely characterise the initial 
field itself. One further reason for pushing measures of £(r) 
to larger and larger scales, is the possibility to correctly 
estimate its zero-point, i.e. the scale on which correlations 
become negative. This is a specific prediction of any vi- 
able model, as it reflects the turnover and peak scale in 
the power spectrum, fingerprint of the horizon size at the 
matter-radiation equivalence epoch. Clearly, as much as 
the weakness of clustering in this regime is a benefit for 
the theory, it represents a hard challenge for the obser- 
vations: statistical fluctuations destroy any possibility to 
detect significant features, such as the zero-point of £(r) , 
if the sample under study does not contain enough objects 
at comparable separations. 

Historically, following the pioneering works of the mid- 
seventies (see Rood 1988 for a review), the industry of 
redshift measurements exploded in the eighties, with the 
completion first of the CfAl (Davis et al. 1982), then of the 
Perseus-Pisces (Giovanelli et al. 1986) and CfA2/SSRS2 
(Geller & Huchra 1989; da Costa et al. 1994) surveys. 
These first large surveys (containing several 10 3 redshifts), 
produced significant advances in the estimation of £ (r) on 
small and intermediate scales (e.g. Davis & Peebles 1983, 
De Lapparent et al. 1988). However, they also showed 
the existence of structures with dimensions comparable 
to their depth. This, while on one side giving rise to spec- 
ulations about the very existence of a transition to ho- 
mogeneity on large scales (see Guzzo 1997 for a recent 
review of this problem) , explicitly indicated the need for 
larger redshift samples. In more recent years, there have 
been a few significant attempts to fulfil this need, as for 
example the surveys based on the IRAS source catalogues 
(see e.g. Strauss 1996), the Stromlo-APM survey (Love- 
day et al. 1992b), and especially the Las Campanas Red- 



shift Survey (LCRS, Shectman et al. 1996). Our contri- 
bution along this direction has been the realization of the 
ESO Slice Project (ESP) galaxy redshift survey, that was 
completed between 1993 and 1996 (Vettolani et al. 1997, 
Vettolani et al. 1998, V98 hereafter). The technical aims 
of the ESP survey were to exploit on one side the avail- 
ability of new deep photometric galaxy catalogues (in our 
case the EDSGC, Heydon-Dumbleton et al. 1989), and on 



the other the multiplexing performances of fibre spectro- 
graphs, as in the specific case of the Optopus fibre coupler 
available at ESO (Avila et al. 1989). The main scientific 
goals were to estimate the galaxy luminosity function over 
a large and homogeneously selected sample with a large 
dynamic range in magnitudes (see Zucca et al. 1997, pa- 
per II hereafter), and to measure the clustering of galax- 
ies over a hopefully fair sample of the Universe. The ESP 
survey, during about 25 nights of observations, produced 
a ~ 85% complete sample of 3342 galaxies with reliable 
redshift. Its combination of depth and angular extension 
is paralleled at present only by the Las Campanas Red- 
shift Survey (LCRS, Shectman et al. 1996), which has the 
advantage of containing a larger number of redshifts. One 
main difference of the LCRS with respect to the ESP, is 
that it is selected in the red, which makes comparison 
with the results presented here particularly interesting. 
One important advantage of the ESP is that it is a purely 
magnitude-limited sample, which makes modelling of the 
selection function easier than in the case of the LCRS, 
where an additional surface-brightness selection was ap- 
plied. 

In this paper we shall discuss the two-point correla- 
tion properties of the ESP survey, both in redshift and 
in real space. Previous or parallel papers, in addition to 
that describing the above mentioned 6j-band luminosity 
function (Zucca et al. 1997), deal with the scaling prop- 
erties (Scaramella et al. 1998), the properties of groups 
(Ramella et al. 1998), and potential biases in the esti- 
mate of galaxy redshifts (Cappi et al. 1998). The survey 
in general is described in Vettolani et al. (1997), while the 
redshift data catalogue is presented in V98. Here we also 
marginally discuss the small-scale galaxy dynamics, for 
its effects on redshift-space clustering. However, a proper 
discussion of the small-scale pairwise velocity dispersion 
and related topics, will be presented in a separate paper 
(Guzzo et al. , in preparation). 

The paper is organized as follows. Section 2 briefly 
summarises the properties of the redshift catalogue. Sec- 
tion 3 discusses the techniques used to estimate two-point 
correlations and their errors. Sections 4 and 5 present the 
redshift- and real-space correlation functions, discussing 
their dependence on luminosity and comparing results with 
other surveys. In section 6 possible biases deriving from 
the observational setup are investigated. The main results 
are summarised in section 7. 
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2. The Survey 

The ESP redshift survey consists of a strip of sky 1° thick 
(declination) by 22° long (right ascension) near the South 
Galactic Pole, plus an additional strip of length 5° situ- 
ated five degrees to the west of the main strip. In total, 
this covers about 25 square degrees at a mean declination 
of -40.25° between right ascensions 22 h 30 m and 01 h 20 m . 
The target galaxies were selected from the Edinburgh- 
Durham Southern Galaxy Catalogue (EDSGC, Heydon- 
Dumbleton et al. 1989), which is complete to bj = 20.5. 
The angular correlation properties of this catalogue were 
studied in detail in Collins ct al. (1992). The limiting mag- 
nitude of the ESP is bj = 19.4, chosen in order to have 
the best match of the number of targets to the number of 
fibres in the field of the multi-object spectrograph Opto- 
pus, mounted at the Cassegrain focus of the ESO 3.6 m 
telescope. More details on the observing strategy and the 
properties of the spectroscopic data are given in Vettolani 
et al. (1997), and in V98. 

The final galaxy redshift catalogue contains 3342 en- 
tries, and is 85% redshift-complete within the survey area. 
At its effective depth (z ~ 0.16), the transverse linear di- 
mensions of the survey are ~ 210 h _1 Mpc by 6.5 h _1 Mpc, 
while its volume is ~ 1.9 • 10 5 h -3 Mpc 3 . 

As can be seen in V98, the survey geometry is fairly 
complex, being composed by two rows of circular Opto- 
pus fields, partly overlapping each other. This results in 
the presence of interstices in which galaxy redshifts were 
not measured, and which obviously have to be considered 
carefully in the analysis of clustering. Also, the complete- 
ness of each Optopus field is not strictly constant (V98). 
This will require a weight to be applied to each object, 
depending on its parent field (see §3.1). The median in- 
ternal error on the redshift measurements is 64kms _1 for 
absorption-line estimates, and 31 kms -1 for emission-line 
estimates (V98). 

In Figure ^ we show a cone diagram of the galaxy dis- 
tribution in the survey. One first qualitative remark to be 
made, in view of the use we shall do of the data in this 
paper, is that the visual impression from this figure sug- 
gests - contrary to more "local" surveys as, e.g. CfA2, a 
typical size of structures which is smaller than the survey 
dimensions. This gives us hope that, possibly clustering 
measures extracted from this survey would be sufficiently 
representative of the general properties of our Universe. 
Clearly, the volume covered by the ESP is still small, due 
to its limited area coverage on the sky, yet the linear sam- 
pling of large-scale structures in two dimensions is un- 
precedented, and matched only by the even larger LCRS. 

Preliminary to the analysis, we performed a series of 
necessary corrections to the raw data catalogue. We first 
applied K-corrections to the observed bj apparent magni- 
tudes in the same way as discussed in detail in Zucca et al. 
(1997). Given the location of the ESP area in the region 
of the South Galactic Pole, we did not apply any correc- 



tion for galactic extinction to the observed magnitudes. 
Observed heliocentric radial velocities were converted to 
the reference frame of the Local Group using the trans- 
formation of Yahil et al. (1977). (We also checked that 
correcting velocities to the CMB reference frame does not 
produce any difference in the results) . We then computed 
comoving, luminosity and angular diameter distances for 
each galaxy, adopting a cosmological model with H Q = 100 
h kms -1 Mpc -1 , and q c = 0.5. 

The correlation analyses are first performed on the 
whole apparent-magnitude limited catalogue. To do this, 
we apply the minimum variance technique (see § 3.1), 
for which we use the survey selection function as com- 
puted in Zucca et al. (1997). We also exclude from this 
analysis those galaxies lying outside of the range of co- 
moving distances 100 < D com < 500h _1 Mpc, to avoid 
using those regions of the sample where respectively ei- 
ther the survey volume or the selection function are dan- 
gerously small. The sample trimmed in this way contains 
2850 galaxies, with a minimum luminosity corresponding 
to Mbj = —15.6 + 51ogh. Using the whole magnitude- 
limited sample is an attempt to extract the maximum sig- 
nal from the available data, but can have some drawbacks, 
as the contribution of galaxies with different luminosities 
is not homogeneous over the sampled scales. 

Therefore, to evidence possible biases and study the 
general behaviour of clustering with luminosity we also 
construct a set of volume-limited subsamples, defined in 
Table 1. Each column lists, respectively (1) absolute mag- 
nitude limit, (2) luminosity distance computed from the 
magnitude limit taking into account the K-correction, (3) 
corresponding comoving distance limit, (4) redshift limit, 
(5) lower cut in comoving distance, (6) effective volume of 
the sample between the two distance boundaries, (7) total 
number of galaxies. Given the small thickness of the ESP 
slice, the lower distance limits are a safeguard to exclude 
that part of the samples where galaxy density is poten- 
tially undersampled by bright galaxies, thus introducing 
shot noise. These low cut-off values have been computed 
as the distance within which, for that absolute magnitude 
limit, one expects less than 10 galaxies within the corre- 
sponding ESP volume, on the basis of the ESP luminosity 
function. The main analyses we shall discuss in the follow- 
ing are based on the first four samples in the table. These 
offer the best compromise between having enough volume 
sampled and a good statistics. However, we also select a 
more "extreme" sample with M < —20.5, which contains 
only 292 galaxies, to study in more detail the existence of 
luminosity segregation at high luminosities. 



3. Computing the Two-Point Correlation 
Function 
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Right Ascension 




-40.8° < 6 < -39.7° 

Fig. 1. Cone diagram showing the large-scale distribution of galaxies in the ESP survey. 
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Table 1. Properties of the volume-limited subsamples analysed. Distances (in Mpc), volumes (in Mpc 3 ) and absolute 
magnitudes (bj system), are computed for h=l, q a = 0.5 
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3.1. Estimators 

There is a long history and literature concerning the the- 
ory of statistical estimators of the two-point correlation 
function £(r). The most widely used one (Davis & Peebles 
1983), is given by 

2n R GG(r) 1 



(1) 



n G GR(r) 

where n G and ur are the mean space densities (or equiva- 
lently the total numbers), of objects in the two catalogues, 
GG(r) represents the number of independent Galaxy-Galaxy 
pairs with separation between r and r + dr, while GR(r) is 
the number of Galaxy-Random pairs, computed with re- 
spect to a random catalogue of points distributed with the 
same redshift selection function and the same geometry of 
the real sample (in other words, a Poisson realization of 
the catalogue). In our case, this implies reproducing the 
exact distribution of Optopus fields on the sky, and apply- 
ing the same selection process as done on the real data. 
One example is the " drilling" of some areas as it was orig- 
inally done on the EDSGC around the position of very 



bright stars ( Heydon-Dumblcton ct al. 1 989j )P|. 

The number of GG(r) or GR(r) pairs in eq. [I] can be 
formed in general including an arbitrary statistical weight, 



as GG{r) = £ ( 



distances of the two objects from the observer (dj and 
dj), and their separation (r = \d\ — dj|). Here the sum 
(i,j) extends over all independent pairs with separations 
between r and r + dr. Any weighting which is indepen- 
dent of the local galaxy density will lead to an unbiased 
estimate of the correlation function, in the sense that the 
mean value of the estimator, if averaged over an ensem- 
ble of similar surveys, would approach the true, under- 
lying galaxy correlation function. However, the variance 
of individual determinations of £ around this mean will 
depend on the particular choice of Wij (r) ; some choices 
lead to more accurate estimates than others. By assuming 
that two-point correlations dominate over higher orders 
(essentially the linear regime) , one finds a useful approxi- 
mation to the minimum variance weights (Saunders et al. 
1992; Hamilton 1993; Fisher et al. 1994a, F94 hereafter), 
as 

1 1 



;(r), that can depend both on the 



i(r) 



(2) 



1 + Aitn G <t>(di)J3{r) 1 + <iim G <l>{dj) J 3 (r) 

where (f>(di) is the survey selection function at the position 
of galaxy i, and Jz(r) = J drr 2 £,(r). The minimum vari- 
ance weights themselves, therefore, require knowledge of £, 
i.e. we have a circular argument for which one must hope 
that an iterative procedure would be convergent. We shall 
instead adopt an a priori model for £ and use it to calcu- 
late the weights; given the host of assumptions employed 
to arrive at (||) , this approach should not be considered as 
too gross an approximation, as also shown, e.g., by F94. 

1 Within the ESP area there are about 10 drill holes with a 
typical diameter of 0.2 deg. 



We model £(r) as a power law, £(r) = (5/r) 1 ' 8 for r < 
30h _1 Mpc, and assume negligible correlations at larger 
separations. This expression results in J 3 (r) = 15.1 r 12 , 
for r < 30 h -1 Mpc, and J3 = 894 at larger separations. 

To test the sensitivity of the results to the adopted 
model, we also compute J3 in a different way, i.e. using 
its definition in terms of the power spectrum P(k). In 
this case, we have Js(r) = (47rr 3 /3) J °° dkk 2 P(k)W t h(kr), 
where Wth(kr) is the Fourier transform of the spherical 
top-hat window function. We compute J3 for the CDM 
power spectrum normalized to unity in spheres of radius 
8h _1 Mpc. [We note that variants of standard CDM, such 
as an open model or a low-density flat model with non- 
zero A, would provide a more realistic power spectrum, but 
as it is clear from the results of the test, this would make 
no real difference for our purposes] . We find no significant 
deviations between the estimates of £(s) computed using 
the two different models for J3, and therefore used the one 
based on £(r) in the rest of the computations. 

The weighting function has to be taken into account 
also when calculating the mean densities n G and ur to be 
used in eq. 0. A full discussion of the merits of different 
estimators for the mean density in a magnitude-limited 
sample is presented in paper II. In general, an estimate of 
the mean density, given a weight for each object and the 
survey selection function, would be of the form 

G f Vs dVw(d)0(d) ' [ ' 

where the sum extends over all galaxies (or random points) , 
while the integral is extended to the whole survey volume 
V s . In practice, we are interested in the ratio of ur and n G , 
so that we only need to estimate for the two sets of points 
the value of the quantity in the numerator (i.e. the effec- 
tive total number of objects). The same arguments leading 
to eq. (|^) result in a similar expression for the minimum 
variance, single-point weights: Wi = (1 + n G cj)(di) J™" 1 ) -1 
(Davis & Huchra 1982); here J™ ax is the maximum value 
of J3. Note that again we have one quantity, n Gl that ide- 
ally should be calculated iteratively. However, since in any 
case our estimator will be unbiased, no matter the choice 
of the weight, we use here - and correspondingly in the pair 
weighting - the mean density as estimated in paper II from 
this same sample, with the same M^, < —12.4 + 51ogh 
cut. 

In our specific case, there is an additional complica- 
tion arising from the varying sampling of the individual 
Optopus fields, as discussed in § |[ This requires the intro- 
duction of another weight, W l , to correct for the missed 
galaxies in each field. Following the same reasoning as de- 
tailed in V98, we define the completeness in each field C l 
as 

NX 



C l = — (4) 

where N l z and N^, are the number of secured galaxy red- 
shifts and the number of objects classified as galaxies in 
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tions: £fe(r), with k — 1, Nb^\ We refer to the set of sepa- 
rations r, , one for each bin of size Ar, as separation space. 
The correlation function £(r) is a vector, £ s in this space 
with components £( r i)- A covariance matrix in separation 
space may be constructed as 

C(n, Tj ) =< (£(r<) - Z(ri)Mrj) - >rea«*, (5) 

where 'Oreaiiz' indicates an average over bootstrap real- 
izations and £ (r^) is the ensemble average of the correla- 
tion function at separation . Since the values of £ at two 
different separations are correlated, C is non-diagonal. 
For this reason one cannot do a straightforward \ 2 fit of 
a model to the observed points. However, C is symmetric 
(C = C) and real, and therefore hermitian, i.e. such that 
it can be diagonalized by a unitary transformation if its 
determinant is non-vanishing. What happens in practice, 
is that the estimated functions are oversampled, so that 
the effective number of degrees of freedom in the data is 
smaller than the number of components of This makes 
C singular, and one has to apply singular value decompo- 
sition (as in F94), to recover its eigenvectors. These define 
the required transformation T, to the orthogonal space, 
i.e. C -> D = T-C-T and £ s -> £ D = T-£ s , where now 
D is diagonal and the components of £ D are independent. 
In this space we can therefore define in the usual way the 
quantity x 2 with respect to our model parameters, and 
minimize it to find their best-fit values. 

4. The Redshift Space Correlation Function 
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the parent photometric catalogue, respectively, in each 
field i, while /* = 0.122 is the observed fraction of misclas- 
sified stars. We therefore multiply the Wi corresponding to 
each galaxy by the weight pertaining to its field, that will 
be given by W l = 1/C 1 . This is not necessary for the 
random sample, that, by construction, is uniform over the 
Optopus fields. 

Eq. (|l|), is based on the definition of £ (e.g. Peebles 
1980), in which it is understood that no represents the 
universal mean galaxy density. Any sample mean density 
will vary about this value with a dispersion at least as 
large as expected from simple Poisson statistics and, in 
fact, clustering on the scale of the survey volume will in- 
crease this dispersion. The precision of the estimate of 
£ at small values (large scales) is limited by this uncer- 
tainty, Suq (the number of points in the random catalog 
iV r = V ■ riR may always be sufficiently increased to elim- 
inate its density as a source of uncertainty, and this is in 
fact what we do here by using typically N r = 100,000). 
Given that ng appears linearly in estimator ([I]) , one would 
conclude that the precision on £ is also linear in Snc- In 
fact, this is not correct, as emphasized by Hamilton (1993) 
and Landy & Szalay (1993), who suggested alternative es- 
timators in which the precision goes as Suq, as would 
seem more appropriate for the second moment of a distri- 
bution. We have computed £(s) using (Q), and with Hamil- 
ton (1993) estimator, and found for our survey equivalent 
results within the errors. This is partly in contrast to what 
was found by Loveday et al. on the Stromlo-APM survey 
(1995), but agrees to the results of Tucker et al. (1997) on 
the larger LCRS (cf. their Figure 1). The results we shall 
show in the following are therefore all based on the Davis 
& Peebles (1983) estimator. 



3.2. Error estimation and model fitting 

Statistical error bars for our estimates of two-point corre- 
lation functions are computed using bootstrap resampling 
as discussed by Ling et al. (1986). F94 have discussed in 
detail the reliability of the bootstrap technique in estimat- 
ing errors for correlation functions; they show that, typi- 
cally, bootstrap errors are a good description of true rms 
errors for separations smaller than ~ 10h _1 Mpc, while on 
larger scales they tend to overestimate them by a factor 
of about 1.5. In the same work, they carefully discuss the 
technique for fitting of models to correlation functions, 
where the data points are not statistically independent. 
We use the same procedure here, which we briefly sum- 
marize in the following. Essentially, we create Nb = 100 
"perturbed" realizations of the sample, by randomly sam- 
pling the original catalog with replacement, (the "boot- 
strapping"). By subjecting each of these realizations to 
the same analysis, we obtain a set of Nb correlation func- 



4-1. Optimal weighting 

To maximise the large-scale signal from the data, we first 
estimate £(s) from the whole ESP redshift catalogue, preparred 
and trimmed as described in § |2| Figure |2| shows the re- 
sult of applying the J3 optimal-weighting estimator to this 
sample, containing 2850 galaxies. Between 3 and 50 h _1 Mpc 
£(s) is well described by a shallow (7 ~ 1.5) power law, 
with a rcdshift-space correlation length s ~ 5h _1 Mpc. 
On larger scales, it smoothly breaks down, crossing the 
zero value between 60 and 80 h _1 Mpc. The bin centered at 
80 h _1 Mpc has indeed a negative value, although it is less 
than ler below zero. On larger scales (r > 100h _1 Mpc), 
there is marginal evidence that the amplitude of £(s) rises 
up again, a behaviour that seems to be shared also by 
other clustering data, as we show below. On scales smaller 
than 3h _1 Mpc, £(s) flattens significantly, so that a sin- 
gle power-law is certainly not a fair description over the 
whole explored range. 

A comprehensive comparison of the best-fit results of 
a power-law £(s) to redshift-spacc correlation functions 
from several surveys is given in Willmer et al. (1998). As 
one could easily guess from our Figure ||, these results 

2 Note that here we shall fit models to one-dimensional quan- 
tities only - i.e. £(s) or projections of £(r p ,7r) 
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Fig. 2. Rcdshift-space correlation function from the whole ESP magnitude-limited catalogue, computed using the 
minimum- variance, ^-weighting scheme. For reference, the dashed line gives an arbitrary power-law function with 
r = 5 h _1 Mpc and 7 = 1.5. 



are strongly dependent on the range of scales over which 
the fit is performed, that often is not the same for the 
different data sets. This makes the comparison of the crude 
best-fit values rather inconclusive. Equivalently difficult 
is any physical interpretation of the measured amplitudes 
and slopes, given that the fit sometimes includes scales 
where power suppression by virialised structures ( "Fingers 
of God") is very effective (s < 2h~ 1 Mpc), while in other 
cases it starts above these, putting more weight on less 
nonlinear scales. Given these ambiguities, here we do not 
perform a formal fit to the ESP £(s) . We shall compute 
a proper fit only to the real-space correlation function, 
for which between 0.4 and 10 h -1 Mpc a power-law shape 
is a rather good description of the data, and where the 
observed shape is at least freed of one major distorting 
effect. 

In Figure |^, we compare the J3 estimate of £(s) from 
the ESP on large scales to those computed with the same 



technique from the LCRS and the Stromlo-APM surveys^ 
The agreement between the three estimates is very good 
between 2 and 30 h -1 Mpc, given also the different galaxy 
selection functions of the corresponding surveys. On larger 
scales we can note that the ESP and Stromlo-APM £(s) 
seem to show slightly more power than the LCRS. This 
evidence for low-amplitude power on large scales, with a 
smooth decay from the power-law shape, resembles that 
observed in the IRAS 1.2Jy rcdshift survey (F94). One 
could probably explain the differences in Figure ^ with the 

3 We limit our explicit comparison to these two surveys, 
as the LCRS is the only other redshift survey with depth 
and angular coverage comparable or superior to the ESP. The 
Stromlo-APM survey, on the other hand, is less deep by ~ 2 
magnitudes and is sparsely sampled, but represents an inter- 
esting comparison because of its large solid angle and the fact 
of being selected from the same photographic material (IllaJ 
plates) . 
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Fig. 3. £(s) on large scales from the ESP, LCRS (Tucker et al. 1997), and Stromlo-APM (Loveday et al. 1992b), redshift 
surveys. The dashed lines show the real-space £(r) as obtained through deprojection of the angular correlation function 
for the APM galaxy catalogue, for two clustering evolution models (Baugh 1996). 



different selection criteria in the three redshift catalogues: 
the LCRS is selected in r (with the additional complica- 
tion of a surface-brightness cut, whose effect is not fully 
clear), and thus should tend to favour earlier types on 
the average, that we know reside preferentially in high- 
density regions. On the other hand, the bj band used for 
the ESP (so as the IRAS infrared band), is expected to 
better select star-forming objects in low-density regions. 
It should probably be expected that ESP galaxies, as well 
as IRAS galaxies, are better tracers of very large-scale, 
low-amplitude fluctuations. 

In Figure ^, we also plot (dashed lines), the real-space 
correlation function obtained by deprojecting the angular 
correlation function w(8) from the APM galaxy catalogue 
(Baugh 1996). This is computed in the case of clustering 
fixed in comoving coordinates (bottom curve), or growing 
according to linear theory (top) . While the details of the 
real-space correlation function will be analysed in the fol- 



lowing section directly from the ESP itself, here we have 
the chance to make already a few interesting remarks. 

First of all, the scale of the break in both £(s) and £(r) 
is consistently indicated by the different surveys to be be- 
tween 50 and 90 h _1 Mpc, with the ESP clearly pointing to 
a value > 50 h -1 Mpc. The agreement between these new 
redshift data and the APM £(r) , directly implies that the 
large-scale power originally seen in the angular correla- 
tion function w(8) from the APM (Maddox et al. 1990) 
and EDSGC (Collins et al. 1992), galaxy catalogues is not 
significantly enhanced by errors in the magnitude scale, as 
claimed by some authors (e.g. Bertin & Dennefeld 1997). 
In fact, Figure || demonstrates that similar power is seen 
in redshift survey data, either based on the same bj plates 
(ESP) or CCD photometry (LCRS). Also concerning the 
detailed shape of galaxy correlations above 3 h _1 Mpc, in- 
dependent data sets show a high degree of unanimity: 
if one ideally extrapolates to larger scales the 'classic' 
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Fig. 4. Redshift— space correlation function from four representative volume-limited samples extracted from the ESP 
as described in Table 1. The dashed line with error bars reproduces the ^-weighted estimate from the whole survey. 
For clarity, points for M < —19 and M < —20 are displayed with a shift in log(s). 



~ — 1.8 slope observed in real space below 3 h _1 Mpc from 
the APM £(r) (but also from previous angular projections, 
see e.g. Davis & Peebles 1983), all surveys are systemat- 
ically above this extrapolation. To reproduce this feature 
(described in the literature as a 'shoulder' or a 'bump', 
see Guzzo et al. 1991, and Peacock 1997), a rather steep 
(~ fc~ 2 ) power spectrum P(k) is required (Branchini et al. 
1994, Peacock 1997). 

A second important point about Figure ^ concerns the 
effect of rcdshift-space distortions. One can see how small 
is the amplitude difference between the redshift-space cor- 
relation functions and the APM real-space £(r) on large 
scales. A linear amplification is indeed expected as a re- 
sult of coherent flows towards large-scale structures. This 
can be approximately expressed as ~ l + ^/3 + i/3 2 , where 
j3 = f^' 6 jb and b is the linear biasing factor (Kaiser 1987). 
Figure || implies, within the errors, a value for this factor 



very close to 1. For example, a value of 1.2, consistent with 
the observed data, would yield (3 = 0.35. 

It is interesting to note that the Stromlo-APM £ (s) is 
the only data set that around ~ 10 h _1 Mpc seems to show 
a significant amplification with respect to £ (r) . (Note that 
this comparison is actually the safest, being the APM cat- 
alogue the parent photometric list of the Stromlo-APM 
redshift survey). The interpretation of this effect is how- 
ever complicated by what we have pointed out in Paper 
II, where we showed how the mean density in this redshift 
survey is about a factor of 2 lower with respect to that 
deduced from deeper samples. The cause for this seems 
to be a negative density fluctuation that we clearly detect 
in the ESP data below z < 0.05, and that is also visible 
in other surveys. This "local" underdensity, is shown to 
be the explanation for both the low normalization of the 
Stromlo-APM luminosity function (Loveday et al. 1992a), 
and the steep number counts observed at bright magni- 
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tudes in several bands (e.g. Maddox et al. 1990). It is not 
trivial to understand how this systematic 50% underesti- 
mate of the mean density would affect the two-point cor- 
relation function, although it clearly reduces by the same 
amount the number of pairs expected at a given separa- 
tion. If this deficit is evenly spread on all scales, then the 
measured £(s) will be practically unaffected, however if 
it corresponds, for example, to a region with larger-than- 
average voids, the net effect will be to boost up £(s) by 
some amount. This would be sufficient to produce the ob- 
served discrepancy between £(s) and £(r) , and further 
demonstrates how extended and accurate the data have 
to be in order to extract dynamical information as the 
value of (3. 

4-2. Volume-limited estimates 

While the ^-weighted estimate has the advantage of max- 
imising the information extracted from the available data, 
it has some important drawbacks that have been not al- 
ways appreciated in previous applications. The main prob- 
lem is that by its own definition it inevitably mixes the 
contribution of galaxies with different luminosities. This 
would not be a problem if galaxy clustering were com- 
pletely independent of luminosity, i.e. if each galaxy traced 
fluctuations independently from its own absolute magni- 



P7 
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tude. As we discussed in §3.1, this technique weighs pairs 
based both on their separation and on their distance from 
the observer. This implies that the main contribution to 
small-scale correlations comes from low-luminosity pairs, 
that are numerous (and dense) in the nearby part of the 
sample and thus better trace clustering at small separa- 
tions. On the contrary, the estimate of £(s) on large scales 
is dominated by the contribution from luminous objects, 
that are in fact the only ones detectable in the distant 
part of the survey. Clearly, if there is any luminosity de- 
pendence of clustering, this technique will tend to modify 
in some way the shape of £(s) . In particular, if luminous 
galaxies are more clustered than faint ones, the J3 weight- 
ing will tend to give a shallower slope for a power-law 
shaped correlation function. Therefore, while this method 
is certainly optimal for maximising the clustering signal 
on large scales (and so our conclusions on the large-scale 
shape of £(s) from the previous section should not be af- 
fected), it might be dangerous to draw from its results 
far-reaching conclusions on the global shape of £(s) , es- 
pecially on small scales. A wise way to counter-check the 
results obtained from the optimal weighting technique, is 
that of estimating £(s) also from volume-limited subsam- 
ples extracted from the survey. In this case, each sample 
includes a narrower range of luminosities and no weighting 
is required, the density of objects being the same every- 
where. This means that in our case Wi(r) = 1, and only 
the W l 's have to be taken into account. 

We have performed this exercise on the ESP, and the 
results for the four main samples defined in Table 1 are 
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Fig. 5. A further test for luminosity dependence of cluster- 
ing in redshift space. £(s) for the brightest (M < —20.5) 
volume-limited subsample of the ESP is compared to one 
of the main samples from the previous figure. Note that 
the scale has been changed to better evidence the differ- 
ence in amplitude on intermediate scales. 



displayed in Fig. || In the figure, we also reproduce the J 3 - 
weighted estimate (dashed line) . To ease visualization, the 
data points for the M < —20 and M < —19 samples have 
been shifted by a constant value (negative and positive re- 
spectively) in log(s). It is clear how for the ESP data the 
J3~based method produces a £(s) which is shallower be- 
low 3 h _1 Mpc. On larger scales the shape is consistent be- 
tween the two methods, and for s > 20h _1 Mpc the opti- 
mal weighting performs better in terms of signal-to-noise, 
than the single volume-limited estimates. The observed 
£(s) is in practice the result of the convolution of the 
real-space correlation function £(r) with the distribution 
function of pairwise velocities. This means that the small- 
scale flatter shape of £(s) in the case of the magnitude- 
limited sample could be in principle produced either by a 
smaller £(r) or by a higher pairwise velocity dispersions. 
We find (Guzzo et al., in preparation), that most of the 
effect is produced by this second factor: the small-scale 
pairwise velocity dispersion between and 1 h" 1 Mpc, is 
higher when measured on the whole survey using the J3- 
weighted method. We measure 1712(1) — 600 kms -1 using 
the J3 estimate from the whole sample, while the volume- 
limited samples give values between 300 and 400 km s -1 , 
depending on the infall model adopted. 

From this plot alone there is no evident sign of a lu- 
minosity dependence of £(s) , since the four estimates are 
virtually the same within the errors. From their analy- 
sis of the SSRS2 sample, Benoist et al. (1996), also find 
that there are negligible signs of luminosity segregation 
for galaxies fainter than M* , which in our case corre- 
sponds to —19.6. They also show, however, that a clear 
effect becomes visible for absolute magnitudes brighter 
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than ~ —20. To check for this effect in more detail, in 
Fig. U we plot £(s) for the brightest sample that can be 
selected from the ESP while keeping a reasonable statis- 
tics, with M < —20.5. This estimate is compared to that 
for the M < —19.5 ~ M* sample from the previous fig- 
ure. Here, we see a significant increase in the amplitude of 
£(s) for the more luminous galaxies, for separations above 
3 — 4h~ 1 Mpc. Even if we do not perform a formal fit for 
the reasons discussed above, we see that now £(s) passes 
through unity around r ~ 8h _1 Mpc, i.e. has an ampli- 
tude very consistent to that observed by Benoist et al. 
(1996) for galaxies with similar luminosity in the SSRS2. 

These authors also point out how the detection of lumi- 
nosity segregation is complicated by volume effects (differ- 
ent samples covering different volumes), and by redshift- 
space distortions. For these reasons, the results of Fig- 
ures H and H are probably meaningful for what concerns 
the shape and amplitude of £(s) on scales > 4 — 5 h _1 Mpc 
(given also the negligible amount of redshift-space am- 
plification that we have shown in Figure ^). On smaller 
scales, however, clustering and dynamics mix up in a non- 
trivial way for the different samples, so that the sim- 
ple redshift-space correlation function is possibly hiding 
a large amount of information. Guzzo et al. (1997) have 
shown how to disentangle some of these effects by studying 
the differences of clustering in real space in the case of the 
Pcrsus-Pisces survey. More recently, Willmer et al. (1998) 
have studied £ (r) and £ (s) in the SSRS2 survey for differ- 
ent luminosity subsamples. We shall see in § || how also in 
the case of the ESP, when studying £(r p ,7r) and its real- 
space projection w p (r p ) , it is possible to evidence some 
weak positive trend of clustering with luminosity also at 
small separations. 



5- tt) and the Real Space Correlation 

Function 

So far, we have considered the correlation among galax- 
ies as located in redshift space, where a map of the galaxy 
distribution is distorted due to the effect of peculiar veloc- 
ities, whose contribution adds to the Universal expansion 
to produce the redshift we actually measure. The effect of 
redshift-space distortions can be shown explicitly through 
the correlation function £(r p ,7r) , where the separation 
vector of a pair of objects is split into two components: 7r 
and r p , respectively parallel and perpendicular to the line 
of sight. Given two objects with observed radial velocities 
v\ = cz\ and v-i = cz 2l following F94 we define the line of 
sight vector 1 = (vi + V2 )/2 and the redshift difference 
vector s = vi — V2 , leading to the definitions 

si > s-s 2 (6) 
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Fig. 6. £,(r p ,ir) f° r the ESP magnitude-limited sample. 
The heavy contour corresponds to £ — 1 ; for higher values 
of £, contours are logarithmically spaced, with A log 10 £ = 
0.1; below £ = 1, they are linearly spaced with A£ = 
0.2 down to £ = 0. The dashed contours represent the 
isotropic correlations expected in the absence of peculiar 
velocities. The map has been Gaussian-smoothed with an 
isotropic filter of width 3 h _1 Mpc. 

Figure |TJ displays the observed £(r p ,Tr) for the ESP 
survey, estimated from the magnitude-limited sample us- 
ing the same minumum weighting technique used for com- 
puting £(s) , but binned linearly into bins of 1 h _1 Mpc. 
As expected, the iso-correlation contours of £(?>,7r) are 
stretched along the line of sight (ir) at small separations, 
with respect to the isotropic case (dashed circular con- 
tours). This is the effect of the relative velocity disper- 
sion of galaxy pairs, <Ji2, and includes the contribution 
of both the large velocity dispersion cluster galaxies, and 
the cooler field population. A careful estimate of <7i2 and a 
discussion of the small-scale dynamics in the ESP survey 
will be presented in a separate paper. Here we limit the 
analysis of £(r p , 7r) to the task of recovering the real-space 
correlation function £(r) . 

The standard way to recover £(r) is to project £(r p , 7r) 
onto the r p axis, in this way integrating out the dilution 
produced by the redshift-space distortion field. The result- 
ing quantity is 



,(7) 



We can then generalize estimator (|l|) to the case of £(r p , ir) 
, now counting the number of pairs in a grid of bins Ar p , 
Att. 



w p (r p ) = 2 dyt{r p ^) = 2 dy £ (r 2 p + y 2 ) 1 ' 2 

Jo Jo 
where the second equality follows from the independence 
of the integral on the redshift-space distortions. In the 
right-hand side of the equation, £ is the real-space cor- 
relation function, evaluated at r = (r^ + y 2 ) 1 / 2 . With a 
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Fig. 7. Results of the fit of a power-law £(r) to the 
projected function w p (r p ) (filled circles), for the whole 
magnitude-limited redshift sample. The inset shows the 
corresponding best-fit values for the correlation length r 
and the slope 7, together with the single-parameter con- 
fidence contours at the 1-, 2- and 3-tr levels. 



power-law model for £(r), i.e. £(r) = (r/ro) 7 the integral 
can be computed analytically, yielding 



w p {r p ) 



r(j)r(V) 



(8) 



where T is the Gamma function. The upper integration 
limit in Eq. (Q), ir up — 00, is chosen in practice so as to 
be large enough to produce a stable estimate of w p . As in 
previous works (Guzzo et al. f997), we find w p (r p ) to be 
quite insensitive to ir up in the range 20 h" 1 Mpc < 7r up < 
25 h -1 Mpc for r p < 10 h" 1 Mpc. In fact, for small values 
of r p , it is important not to choose a too small upper limit, 
since this would leave out small-scale power that is present 
at large 7r's due to the stretching. On the other hand, the 
choice of too large a value for tt up , has the only effect of 
adding noise into w p . 

The computed w p (r p ) for the whole magnitude-limited 
ESP catalogue is shown in Figure [?] (solid circles), together 
with the results of the fit of a power-law model for £(r) 
(solid line and inset). Error bars are as usual computed 
over 100 bootstrap realizations of the sample. The corre- 
sponding best-fit values are r a = ^.lh^^h^ 1 Mpc and 
7 = 1.67^0 09- These values are smaller than those mea- 
sured from the LCRS, for which Jing & Borner (1996) 
obtain r Q = 5.06 ±0.12 and 7 = 1.862 ±0.034. The reason 
for this difference is probably to be found in the different 
galaxy populations that these two surveys preferentially 



select, an effect possibly enhanced by the faint limiting 
magnitude of the two surveys. In fact, from the blue- 
selected SSRS2 (limited to to = 15.5, i.e. 4 magnitudes 
brighter than the ESP), Willmer et al. (1998) measure a 
real-space correlation function with r a = 5.36 ± 0.33 and 
7 = 1.86 ± 0.16, i.e. virtually identical to the red-selected 
LCRS. Clearly, given its fainter magnitude limit, in the 
ESP a much larger population of intrinsically faint, possi- 
bly less clustered, galaxies enters the sample. This could 
well be the reason for the smaller global clustering ampli- 
tude measured from the magnitude-limited sample, where, 
as we explained above, faint nearby galaxies play a domi- 
nant role on small scales in the J3-weighted technique. 

Let us now look at £(r p , w) for the four "best" volume- 
limited samples, whose plots are shown in Figure b[ The 
four panels display a number of interesting features. First 
of all, one immediately sees the striking compression of the 
contours at large (> 10h _1 Mpc) r p s for the two fainter, 
smaller-volume samples in the top panels. This effect is 
similar to what one would expect as consequence of infall 
motions onto superclusters, and if real could be used to 
estimate again the parameter (3 that we discussed in the 
previous section (Kaiser 1987). A similar effect is not seen, 
however, in the two more luminous, deeper samples. An 
explanation for this behaviour can be attempted by look- 
ing at the cone diagram of Figure [j], and considering the 
limiting distances for the different volume-limited sam- 
ples, (Table 1). In the case of the —18.5 and —19.0 sam- 
ples (z < 0.11 and z < 0.13, respectively), one can see that 
two superclusters, nearly perpendicular to the line of sight, 
are dominating the galaxy distribution. This is clearly en- 
hanced by their position near the far end of the sample, 
where the sampled volume is maximum. In other words, 
more nearby, differently oriented, structures are just "cut- 
through" by the thin ESP slice, and weigh less, while these 
two are included nearly in their full size. In this way they 
contribute an excess number of pairs along r p , and thus an 
anisotropy (compression) of the contours, which has noth- 
ing to do with a true infall. Note indeed how the effect is 
maximum for the M < — 19 sample, when both struc- 
tures are fully included within the sample limit. Despite 
the depth of the ESP, therefore, these two subsamples do 
not yet provide an unbiased (with respect to direction) set 
of superstructures. Apparently, a more isotropic situation 
is reached within the two brighter samples, that cover a 
larger volume. In these two cases, however, the signature 
of large-scale motions is very weak. An accurate study of 
this kind of distortions would require a careful modelling 
of all selection effects, and even with this a measure of f} 
from the ESP looks rather difficult (see e.g. Fisher et al. 
1994b for details on the technique). 

Still looking at Figure ||, one can also notice that the 
small-scale stretching of the contours is more pronounced 
for the more luminous galaxies (M < —20), than for 
the fainter ones. Again, volume effects mix together with 
possible luminosity effects, but one interpretation of this 
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Fig. 8. £(r p ,7r) for the four main volume-limited subsamples of the ESP catalogue. Contour levels are defined as in 
Figure |[ 



could be that luminous galaxies are more likely to be found 
in dense, high-velocity dispersion regions. In fact, this is 
directly confirmed by our analysis of the luminosity func- 
tion of galaxies in ESP groups and clusters (Ramella et al. 
1999), where we find that galaxies in high-density regions 
show a brighter M* with respect to "field" objects. It is 
therefore natural that selecting luminous galaxies, we also 
select galaxies with higher pairwise dispersion. 



This picture is perfectly consistent with the behaviour 
of the real-space correlation function for these samples. 
Figure ^| shows the results of the fit of a power-law £(r) 
to the projected function w p (r p ) , while Table |2| sum- 
marises the corresponding best-fit values and errors. In- 
terestingly, in real space there is a weak, but systematic 
trend in both r Q and 7 towards larger values for larger lu- 
minosities. This is consistent with the higher small-scale 
velocity dispersion inferred from the stretching of the con- 
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Fig. 9. Best-fit results of a power-law £(r) to the projected function w p (r p ) for the same four volume -limited samples 
of the previous figure. Left panels give the observed w p (r p ) with the corresponding best-fit power law. Right panels 
show the best-fit r Q and 7 values, together with their 1-, 2-, and 3-cr confidence levels. 



tours of £,(r p ,ir) . One can also note from the left panels 
of Figure ^, that the effect is mostly produced on scales 
smaller than 2 h^ 1 Mpc, while no big change in amplitude 
is evident at larger separations. Here we understand why 
no evident trend could be seen for the same four samples 
when we compared their redshift-space clustering through 
£(s) : more luminous galaxies are indeed more clustered 
on small scales [higher £(r) ], but their pairwise dispersion 



increases accordingly and so the damping of £(s) . The net 
result is that in redshift space the two effects nearly cancel 
each other, and indeed nearly no sign of luminosity segre- 
gation is typically seen up to luminosities corresponding 
to M ~ -20 (Benoist et al. 1996). 

Similar real-space clustering analyses have been per- 
formed in the past on more "local" samples as the Perseus- 
Pisces (Guzzo et al. 1997), and the SSRS2 (Willmer et al. 
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1998) redshift surveys. Even allowing for systematic differ- 
ences in the Zwicky and bj magnitude systems, in general 
the ESP estimates tend to be smaller than those from 
samples of comparable mean luminosity extracted from 
these two surveys. In particular, values of the correlation 
length measured here for generic galaxies are comparable 
to those of spiral-only samples of similar luminosity from 
the Perseus-Pisces survey. This would again suggest that 
the ESP is on the average richer in late-type galaxies. 



6. Tests for Potential Pathologies in the ESP 

The technical characteristics of the Optopus fibre spec- 
trograph used for the redshift survey, introduced a bias at 
small angular separations, due to the finite size of the fibre 
connectors. This forced a minimum separation of 25 arc- 
sec, below which two objects could not be observed with 
the same plugged plate. Galaxies which are closer on the 
sky, therefore, had to be observed in two different expo- 
sures of the same Optopus field, to be able to collect a 
spectrum for both of them. For some of the survey fields 
(the most crowded) , these repeated observations were not 
possible. The consequence, explicitly shown in V98f], is 
that when plotting the distribution of separations on the 
sky between each galaxy and its nearest neighbour, this is 
significantly more skewed towards small separations for 
the objects without redshift, than for the total photo- 
metric parent sample. This could in principle introduce a 
bias into the correlation function estimate, and the conse- 
quences for £(s) and £(r p , tt) have therefore to be explored 
carefully. 

First, let us simply consider the fact that at the peak 
of the global redshift distribution function [oc z 2 <j){z)\, 
around z ~ 0.1, the minimum fibre separation corresponds 
to less than ~ 0.1 h _1 Mpc. This should in principle affect 
directly only the first bin we used in our estimates of cor- 
relation functions: for r much larger than the pair separa- 
tion, i.e. the large scale bins, the proper pair count statis- 
tics should be guaranteed by the area-weighting scheme. 

To test any residual effect directly, we have made a 
further, simple exercise. Undoubtedly, the maximum bias 
on the correlation function occurs if all the missing ob- 
jects do really lie at the same distance of their nearest 
neighbour. Let us therefore make the extreme assumption 
that all the close angular pairs with one missing redshift 
are true physical pairs, i.e. that cz\ ~ cz<i. This situa- 
tion should at least bracket the maximum bias we can 
expect from the missing pairs. To reproduce this case in a 
realistic fashion, we have selected all objects without red- 
shift that are closer than 50 arcsec to a companion with 



4 Actually, Figure 3 of that paper shows that the bias is 
present out to separations of ~ 50 arcsec, indicating that in 
addition to the mechanical limitation, there has also been some 
systematic tendency to avoid the very close pairs, when posi- 
tioning the fibres during the observations. 
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Fig. 10. Test for the fibre distance bias. Open circles show 
the estimate of £(s) from the modified sample where all 
close companions without redshift are assigned the red- 
shift of the nearest neighbour (see text), compared to the 
estimate of £(s) from the original catalogue. Both esti- 
mates use the J3 optimal weighting method. 



measured redshift cz\. We have then assigned to each of 
these objects a fake redshift CZ2 = cz\ + S v , where S v 
is a displacement extracted at random from a Gaussian 
distribution with ax% = 350kms _1 0. Through this pro- 
cedure, the resulting sample with "measured" redshifts 
contains about 100 additional galaxies. On this sample we 
repeated our computation of £(s) and £(r p ,7r) ■ The J3- 
weighted redshift -space correlation function is compared 
in Figure [To] to that from the original sample. Evidently, 
the addition of the nearby pairs has very little effect, indi- 
cating that the minimum-separation bias in the fibre po- 
sitions has no consequences on our correlation estimates. 
Analogously, we re-computed w p (r p ) and obtained virtu- 
ally no change in the best fitting parameters for £(r) . 



7. Summary 

The main results obtained in this work, can be summarised 
as follows. 

1. The redshift-space correlation function £(s) for the 
whole magnitude-limited ESP survey can be grossly 
described by a power law with 7 ~ 1.5 between 3 
and ~ 50h _1 Mpc, where it smoothly breaks down, 



5 This value for CT12 can be taken as a fiducial value for the 
pairwise dispersion in the field, as shown by Guzzo et al. (1997); 
note that in reality the line-of-sight distribution of pairwise 
velocities is best fitted by an exponential function, rather than 
a Gaussian. This would however simply increase the tails of the 
simulated displacements: by using a Gaussian we are placing 
the pairs a bit closer on the average, and therefore testing the 
effect on £(r) in an even more strict way. 
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Table 2. Summary of the best-fit parameters of the real-space correlation function obtained from w p (r p ) for the 
different samples analysed. 



crossing the zero value between 60 and 80 h -1 Mpc. 
However, if one considers the overall shape from 0.1 to 
100 h _1 Mpc, a power law is a poor description of the 
observed £(s) , and it is not very meaningful to arbi- 
trarily select one range or the other for a formal fit. 
In particular, below 2h _1 Mpc, £(s) has a much shal- 
lower slope (7 < 1.) . This is shown to be mostly due 
to redshift-space depression by virialized structures, 
whose effect is enhanced by the J3 optimal-weighting 
technique. 

Comparison to the only other redshift survey which 
provides a combination of depth and angular cover- 
age similar or superior to ESP, i.e. the LCRS, shows 
a very good agreement above 2 h _1 Mpc. Minor differ- 
ences are most easily explained in terms of the different 
selection procedures of the two galaxy samples. Also 
comparison with a survey with very different geometry, 
depth and redshift sampling (but with the same magni- 
tude system), the Stromlo-APM, shows a very similar 
behaviour for £(s) . The unanimity of these surveys 
with rather different geometries and selection criteria, 
shows that the shape and amplitude of clustering for 
~ M* optically-selected galaxies is now rather confi- 
dently established, at least between 1 and 20h _1 Mpc. 
Comparison of these redshift-space correlation func- 
tions to the best available real-space correlation func- 
tion from the APM survey, shows a negligible ampli- 
fication from streaming motions, implying a low value 
for the parameter j3 = Q®- 6 /b. 

Computing £(s) from volume-limited subsamples of 
the ESP, we find no significant luminosity segregation 
in redshift space up to limiting absolute magnitudes 
as bright as Mbj = —20. We start seeing some effect 
only for a subsample with M < —20.5, for which the 
amplitude of £(s) increases by a factor ~ 2 on inter- 
mediate scales. This global behaviour in redshift space 
is in general agreement with what found by Benoist et 
al. (1996) on the SSRS2. 

These results also evidence how the small-scale shape 
of £(s) as estimated using the ^-weighting technique 
is systematically flatter than that observed from any 
of the volume-limited samples, where no weighting is 
used. Apparently, this is produced by enhanced weight 
put by the technique on pairs of galaxies with high rel- 
ative velocity (i.e. within nearby clusters), resulting in 



a high pairwise velocity dispersion and consequently in 
a stronger damping of £(s) below 2 h _1 Mpc. Although 
the effect is observed to be more severe in the case of 
the ESP sample, probably because of its thin-slice ge- 
ometry [the LCRS £(s) below 2 h _1 Mpc is closer to the 
volume-limited estimates], these differences underline 
the importance of comparing weighted and unweighted 
estimates. 

5. Studying clustering in real space through £(r p , n) and 
its projection w p (r p ) , we can evidence a mild de- 
pendence on luminosity also for magnitudes fainter 
than —20.5. This effect is mostly confined to separa- 
tions r < 2h~ 1 Mpc. The shape of the contours of 
£(r p ,7r) also clearly shows how the small-scale pair- 
wise velocity dispersion, and the corresponding distor- 
tions (mostly produced within virialised structures), 
are increased in more luminous samples: more lumi- 
nous galaxies inhabit preferentially high-density re- 
gions, where they consequently have to move faster 
with respect to each other. This is the specific reason 
why no evident dependence on luminosity is observed 
on small scales in redshift space, where the two effects 
cancel each other. One could speculate that this mild 
effect evidenced on small scales might be the product 
of dynamical evolution of galaxies within virialized re- 
gions, while that observed only for highly luminous 
galaxies on large scales, given the time scales implied, 
can only be produced ab initio in a biased process of 
galaxy formation. 

6. Fitting w p (r p ) for the whole magnitude-limited ESP 
catalogue with a power-law £(r) , we obtain as best 
parameters r Q = 4.15j:°J°/i" 1 Mpc and 7 = 1.67±g;g^. 
These amplitude and slope are slightly smaller than 
those from most other optical surveys, suggesting that 
late-type galaxies are probably slightly favoured by 
the ESP selection, with respect to early types. 
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